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Definitions  /  Observations 
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pcf  -  P(t  <tL)=  Fj  (tL  )  Prob.  of  Time  to  F ailure 
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Definitions  /  Observations 


(“Hard”  Failures)  (“Soft”  Failures) 


Quality:  Performance  Reliability  at  t  =  0 

OR 

Conformance  to  Specifications  at  t  =  0 


Performance  Reliability  =  Quality  over  Time 
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Design  for  Lifecycle  Cost 
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C,  (d,  X,  I  ,  r)  =  C,(d,  x)+  C,  (d,  X,  t, ) + C,f  (d,  X,  t  r) 

7  /  \  \ 

Lifecycle  Production  Inspection  '  Expected 

Cost  Cost  Cost  Variable  Cost 


Final  time  Interest  rate 

,  >  I 

Cy(d,'K,tf,r)=  \  c F(t)e~rt f j  {t)dt 

°f  \ 

Cost  of  failure  PDF  of  time  to 

at  time  t  failure 

Fj  (tL  )  =  P(3t  e  [O,  tL  ],  such  that  g('X(t),t)  <  0) 
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Design  for  Lifecycle  Cost 
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min  CL  (d,  X,  tf ,  r) 

d,jiX’°x 


#  of  time-dependent 
limit  states 


s.  t. 


n 


Fj  {d,X,t)<  p\  for  i  =  1,2,- 

l  J  l 
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e.g.  p  f  -  0.05  after  t  =  5  years 
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SPL  (tj) 


Calculation  of  Cumulative 
Probability  of  Failure 


■  Series  system  approach  for  [0,  tf  ] 

■  Reliability  estimation  at  time  tL  ( 0  <  tL  <  tf ); 
multiple  MPP  case 
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Calculation  of  Cumulative 
Probability  of  Failure 


Fj  (tF )  =  P(3f  e  [(),  tF  \  such  that  g(x(t),t)<  0) 


At 


0  t\  ^2 


t 
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-2  ^K- 1 


t, 


tF  =  KAt 


time 
- ► 


Series  System 
Prob.  of  Failure 
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Cumulative  Probability  of  F ailure 

(PHI2  Method) 


up-crossing  rate 
r  i  f 

£[AT(0,^)J=  \v+(t)dt 

/ 

#  of  up-crossings  in  [0,  tF  ] 
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(?)=  lim 


p[(g(x(t), t)>  0)n (g(x(f  +  At), t  +  At)<  0)] 
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Cumulative  Probability  of  Failure 
(Set  Theory  Approach) 

time 
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Cumulative  Probability  of  Failure 
(Set  Theory  Approach) 
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/>(AF, ) «  P(CF,  u  CFW )- />(CF, ) 


Probabilistic  Re-Analysis  MC  Method 
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Cumulative  Probability  of  F ailure 

(PHI2  Method) 
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Xl  ~7v(2.29*10“3,0.229*10"3) 

X2  ~  iv(2*1011,0.2*1011) 


X3  ~  iV(0,0.7) 
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Cumulative  Probability  of  Failure 

(PHI2  Method) 
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Correlation  Coefficient, 


Cumulative  Probability  of  F ailure 

(PHI2  Method) 
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Cumulative  Probability  of  F ailure 

(PHI2  Method) 
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SPL  (tj) 


Reliability  Estimation  at  Time  tL 


> 
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■  Time-dependent  reliability 

■  Multiple  MPPs 

■  Niching  genetic  algorithm 
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Main  Points 


>  Time-dependent  reliability  problem  is 

transformed  into  a  time  independent  reliability 
problem  for  solving  a  level-crossing  problem 


>  Multiple  disjoint  failure  domains  may  exist 
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Time-dependent  Reliability 
Analysis  of  a  Random  Process 


>  Transfer  a  time-dependent  process  to  a  time  independent  response: 

S  =  S(X,t),  tirin<t<t^  =>  5imx(X)  =  max  S(X,t) 

t  •  <  t<t 
min  ~  max 

>  Time-independent  limit  state: 

g(x)  =  5allowable-5max(X)  =  0 
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Time-dependent  Reliability 
Analysis  of  a  Random  Process 

>  Double  loop  approach: 

Outerloop :  Find  MPP' s 
P,  =  min|U|| 

s.t.  g(x)  =  sattowabk(X)  -  5raax  (X)  =  0 

X  =  [ix  +  oxU 

Inner  loop : 

Find  maximum  response 
in  time  domain 
S™x(X)=  max  SAX,t) 

J  ^min  —  t  —  ^max 

End  of  innerloop 
End  of  outerloop 
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Time-dependent  Reliability 
Analysis  of  a  Random  Process 


Oakland 
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>  Definitions 

S  Actual  MPPs 
S  Niches 

>  Observations 

S  Niche  center  is  an 
approximate  MPP 

S  Niching  GA  finds  ALL 
approximate  MPPs 
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Time-dependent  Reliability 
Analysis  of  a  Random  Process 

>  Identification  of  Actual  MPPs 

■S  Linearized  limit  states 

G'„  (U)  =  (u  - MPP  )■  MPP  =  0 


or 


gpl(u)=(u-mppJ-vg„  =0 

S  Collinearity  criterion  for  actual 
MPP 

/  X  MPP  VG 

cos  Z(MPPp ,  VGp )  = - - p—  =  ±1 
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MPP  VG 

r  t 

S  Practical  criterion 

>i  —  p 
—  1  6  MPP 


MPP/; 

VGP 

MPP/; 

VGP 

0  ^MPP  ^  1 


G,(U)  = 
Safe  domain 
G(U)  >  0 
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Time-dependent  Reliability 
Analysis  of  a  Random  Process 
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>  Identification  of 
Independent  and 
Significant  MPPs 

>  Two  MPPs  are 
independent  if  the 
correlation  coefficient 
is  very  close  to  +1 

>  Reliability  index 
identifies  significant 
MPPs 


Ppq  =  cosZ(MPPp,MPPj 


MPP/? 

MPPp 


MPP 

H 

MPP 
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^8  Time-dependent  Reliability  Analysis  of  a 

Random  Process 


Oakland 


UN\fR5TTY 


>  Estimation  of  Probability  of  Failure 

S  Convex  safe  domain  formed  by  linearized  limit  states 


M  p  : 

Gp  (U)  =  U  •  MPP);  -  MPP„  •  MPP„ 
<0  for  p  =  1,2,...,  M 

S  Second  order  (bi-modal)  bounds 

M  f  p~ 1  A 

pft  +Emax 

P= 2 


#  = 


MPP. 


(  P~X  (  \  1  MM 

pP ~ZpiMP 0  * Pf  ^ X Pf,  -Zmax p\mp nMJ 

V  <7=1  J  P= 1  P= 2  ?<p 


^  MCS  on  the  linearized  convex  safe  domain 


MM 

O'  =|J;V/ij  =|J(u-MPPp  MPP  MPP(,  <0) 
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A  Niching  GA  Method  for  Identifying 

Multiple  MPPs 

>  Greedy  Fitness  Sharing 

S  Fitness  sharing  is  a  popular  niching  technique. 

■  Decreases  chance  of  mating  in  densely  populated 
niches 

■  Increases  chance  of  mating  in  sparse  niches. 
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A  Niching  GA  Method  for  Identifying 

Multiple  MPPs 
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S  Fitness  function  for  identification  of  MPPs  using  RIA 
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Numerical  Example 
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•Model  Description 

>3%  critical  modal  damping 
>  Carried  out  for  0.1  seconds. 


P(t)  =  25  cos(300;rw(/  -  0.004))  Ibf 


pit)  =  cos(500/tf)  psi 


max 
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Numerical  Example 

Plate  properties 


Mean 

Mx 

Cov  = 

L1 x 

Std.  Variation 

Length  (a) 

2in 

0 

0 

Height  (b) 

2in 

0 

0 

Thickness 

0.0877in 

0.0513 

0.0045in 

Weight  Density 

0.282 

lbs/in3 

0 

0 

Mass/Weight  Factor 

2.59  X10-3  sec2/in 

0 

0 

Young’s  Modulus 

30.0  X106  lbs/in2 

0.05 

15  X  10s  lbs/in2 

Poisson’s  Ratio 

0.3 

0 

0 
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Numerical  Example 
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w 
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CM 


U1  "  t 
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Multimodal  maximum 
displacement  response 


c  max 

S  =  w 


max 


MPP 


4* . 


yG(U)  >  0 
g3(U)  =  0  ^  mpp< 


MPP 


I.*.-'' 


G(U)  <  0 


.  O 

MPP?  G2(U)  =  0 


G(U)  >  0  f 
G1(U)  =  0 


G(U)  <  0 


0.3 

0.2 

0.1 

-10 


—  Limit  State  : 

G(U)=HmLS -0.3  =  0 


10 
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Numerical  Example 
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Identification  of  independent  and  significant  MPPs. 

>  Form  a  linear  approximation  of  each  limit  state  at  each 
MPP. 
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Numerical  Example 
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Estimation  of  probability  of  failure  using  the  Ditlevson’s 
second-order  bounds: 


pLf  =  puf  =  3.3298  xlO"4 


MC  simulation  with  1,000,000  samples  using  linearized 
safe  domain  yields: 

PfLMCS=3.51xlO-4 


If  MPP4  is  removed:  pf  =  3.3287x10  4  i=> 
i==>  MPP4  is  insignificant 
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Observations  /  Future  Work 

■  Calculation  of  cumulative  prob.  of  failure 
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>  Series  system  approach  for  [0,  tf  ] 

>  Reliability  estimation  at  time  tL  ( 0  <  tL  <  tf ) ; 
multiple  MPP  case 


Response  (t) 


Deterministic  Re-Analysis  method 
using  PROM,  MCA  and  Kriging 

Time-dependent  metamodeling 


■  Series  system 
probabilistic  analysis 


Probabilistic  Re- Analysis  methot 
using  single  MC  simulation 


■  Efficient  solution  of  design  optimization  problem  for 
lifecycle  cost 
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Q&A 


**Disclaimer:  Reference  herein  to  any  specific  commercial  company,  product,  process,  or  service  by  trade  name,  trademark, 
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United  States  Government  or  the  Department  of  the  Army  (DoA).  The  opinions  of  the  authors  expressed  herein  do  not 
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endorsement  purposes.** 
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